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The heavy quark collisional scattering on partons of the quark gluon plasma (QGP) is studied in a QCD 
medium at finite temperature and chemical potential. We evaluate the effects of finite parton masses and widths, 
finite temperature T and quark chemical potential fiq on the different elastic cross sections for dynamical quasi¬ 
particles (on- and off-shell particles in the QGP medium as described by the dynamical quasi-particles model 
“DQPM”) using the leading order Born diagrams. Our results show clearly the decrease of the qQ and gQ 
total elastic cross sections when the temperature and the quark chemical potential increase. These effects are 
amplified for finite fiq at temperatures lower than the corresponding critical temperature Tc{liq). Using these 
cross sections we, furthermore, estimate the energy loss and longitudinal and transverse momentum transfers of 
a heavy quark propagating in a finite temperature and chemical potential medium. Accordingly, we have shown 
that the transport properties of heavy quarks are sensitive to the temperature and chemical potential variations. 
Our results provide some basic ingredients for the study of charm physics in heavy-ion collisions at Beam 
Energy Scan (BES) at RHIC and CBM experiment at FAIR. 

PACS numbers: 24.10.Jv, 02.70.Ns, 12.38.Mh, 24.85.+p 


I. INTRODUCTION 

The exploration of the phase diagram of strongly interacting matter is a major field of modem high-energy physics. The tran¬ 
sition from hadronic to partonic degrees of freedom at high temperatures or high baryon densities is one of the most interesting 
challenges of relativistic heavy-ion physics. The discovery of this phase transition is expected to elucidate some of the funda¬ 
mental aspects of Quantum Chromo Dynamics (QCD), i.e. confinement and chiral symmetry breaking. Moreover, the different 
phases of the QCD phase diagram play an important role in the evolution of the early universe and the structure of the core of 
neutron stars [1]. 

The transition from a hadronic to a partonic medium at small net-baryon densities is known to be a crossover. On the other 
hand, at high baryon densities one expects new phases of strongly interacting matter [2] . In particular, a first-order deconfinement 
phase transition with a critical endpoint or a chiral symmetry restoration without deconfinement - leading to a quarkyonic phase 
- may occur at larger baryon densities. 

Recent results from a statistical analysis of particle ratios measured in Pb+Pb and Au+Au collisions at Schwerlonen- 
Synchrotron (SIS), Alternating Gradient Synchrotron (AGS), Super Proton Synchrotron (SPS) and Relativistic Heavy-Ion 
Collider (RHIC) energies point in this direction [3-7]. The phase boundary between quark-gluon and hadronic matter and 
the location of the (possible) critical endpoint is suggested by lattice QCD calculations [8, 9] to occur for values of {Jb larger than 
about 400 MeV; for jis smaller than 400 MeV one expects a smooth cross over from the hadronic to the partonic phase. 

Besides this hypothetical first order transition a rich phase structure might occur at high baryon chemical potential {Jb or quark 
chemical potential {Jq. 

Previous studies have shown that the high density regime of QCD is accessible with heavy-ion collisions at moderate collisions 
energies [10]. This is supported by the study of the freeze-out points from hadron-gas models as a function of the temperature 
and chemical potential [II] or hadronic transport models like the hadron-string-dynamics (HSD) approach from [12, 13], where 
the highest baryon densities are predicted for moderate collisions energies. 

Recent progress in the exciting field of QCD at high baryon densities is driven by new experimental data from the Beam 
Energy Scan (BES) program at RHIC (y(s^=7-200 GeV) and NA61 at the SPS (^/s^=6A-llA GeV). In future experiments 
of the Compressed Baryonic Matter (CBM) collaboration at PAIR (-y/^a^a^= 2.7-8.3 GeV) and the Multi-Purpose Detector (MPD) 
at NIC A (y(s^=4-II GeV) will provide additional information. The aim of all these experiments is to explore the QCD phase 
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diagram at high net baryon densities and moderate temperatures in nucleus-nucleus collisions. This approach is complementary to 
the studies of matter at high temperatures and low net baryon densities performed at RHIC and the Large Hadron Collider (LHC) 
which are designed to study the properties of the deconfined medium at the highest available energy densities (and temperatures). 
Since lattice QCD (IQCD) does not provide robust results in this regime one has to rely on effective QCD models that match 
IQCD for \iq — 0. Furthermore, one has to look at the characteristics of the medium systematically in terms of collision energy 
and system size, by studying the strangeness, charm, collective flow and fluctuations possibly on an event-by-event basis. 

Charm physics is one of the promising signals that can be studied at upper FAIR energies. Indeed, both hidden and open charm 
are expected to contribute to the total charm production also in this energy range. The total charm cross section experimentally 
is not well known close to threshold and the predictions for A-fA rely on parametrisations of experimental data and imply large 
uncertainties towards threshold. Note that also perturbative QCD (pQCD) calculations show large uncertainties [14]. Therefore, 
close to the kinematic production thresholds, an unknown territory has to be systematically explored [15, 16]. Furthermore, 
experimental data close to threshold for the elementary charm cross sections have to be taken before more robust conclusions can 
be expected for p-^A and A-fA collisions at FAIR/NIC A energies. Note that there are presently no p-^N data below 20 GeV and 
no A-fA data below top SPS energies for charmonia. 

The anomalous suppression of charmonium due to screening effects in the Quark-Gluon Plasma (QGP) was predicted to be 
an experimental signal of the QGP formation by Matsui and Satz [17] since particles containing heavy quarks like charm are 
produced predominantly in the early stage of the collision due to high kinematical thresholds. Heavy-flavour physics at FAIR 
thus aims to explore how charm is produced at beam energies close to the kinematical thresholds and how charm propagates in 
hot and cold nuclear matter as well as in the partonic phase. Besides, the production mechanisms of D and // i/r mesons will be 
sensitive to the conditions inside the early fireball. 

The scenarios for charm production in A-fA collisions based on either hadronic [13, 18] or partonic (statistical hadronization 
model, SHM) [19] models have given different predictions. Especially the ratio of hidden to open charm (//i/r/D) from the 
hadronic HSD model [15, 16] differs substantially from the one in the partonic SHM model [19] since the //i/r/D ratio depends 
on the energy in the c.m.s. (y^) in the first model and is yT independent in the second one. Besides, the JjD ratio is about 
one order of magnitude higher in the hadronic than in partonic production scenarios which is due to the much lower threshold 
for //i/r production in N-fN collisions than for D + 5 pairs. Therefore, the charm production is sensitive to the phases of matter 
and the ratio of hidden to open charm appears as a very promising probe for the production of charm and its propagation in the 
medium. In addition the collective flow of charm is expected to provide valuable information on the interaction strength of charm 
with its medium being of hadronic or partonic nature. Furthermore, it is generally questioned that charm degrees of freedom 
might achieve a chemical equilibrium in the fireball such that microscopic transport approaches are mandatory to shed some light 
on the nonlinear charm dynamics. 

With the future aim to implement the charm dynamics non-perturbatively into the Parton-Hadron-String Dynamics (PHSD) 
approach we have to specify the interaction cross sections of the charm degrees of freedom (g) with the light partonic degrees of 
freedom incorporated in PHSD, i.e. dressed quarks, antiquarks {q^q) and gluons (g). Accordingly, in this study we will compute 
the off-shell cross sections for the reactions qQ qQ and gQ gQ taking into account the quasi particle nature of the quarks 
and gluons at finite T and flq in PHSD which are adopted from the dynamical quasi-particle model (DQPM) [20-22] . These cross 
sections are then used to evaluate the energy and momentum losses of a heavy quark as a function of heavy-quark momentum, 
temperature and quark chemical potential. Our results Anally will provide the basic ingredients for the microscopic study of 
charm physics for the Beam Energy Scan (BBS) program at RHIC and the future CBM experiment at FAIR. For a first step in this 
direction we refer the reader to Ref. [23] where some of the quantities have already been evaluated at = 0 at finite temperatures 
relevant at RHIC or LHC energies. 

The paper is organized as follows: We first present in section II the basic ingredients needed for the heavy quark scattering in 
a finite temperature T and quark chemical potential jlq medium. Therefore, we fix the coupling constant, the parton masses and 
spectral functions and the gluon and fermion propagators as given by the DQPM at finite T and jlq. The analysis of the on- and 
off-shell kinematics and the calculations of the on- and off-shell elastic cross sections of the scattering of a heavy quark (g) and 
light quark {q) and gluon (g) in a partonic medium at finite T and jlq are specified in Section III. In Section IV, we calculate the 
interaction rate of the heavy quark in a such medium. Furthermore, we will perform a quantitative analysis of the heavy-quark 
energy loss (sec.V) and momentum loss (sec. VI). Throughout Sections IV-VI we will point out the effects of finite masses and 
width, finite temperature and chemical potential on the heavy-quark transport properties. In Section VII we summarize the main 
results and point out their future applications. 


II. DYNAMICAL QUASI-PARTICLES AT FINITE TEMPERATURE AND CHEMICAL POTENTIAL 

The scattering of heavy quarks in vacuum and in a QGP medium in lowest-order QCD perturbation theory (pQCD) has 
extensively been studied in the literature [24, 25]. Recent developments reconsidered the concept of perturbatively interacting 
massless quarks and gluons as constituents of the QGP, which scatter according to the leading (Bom) diagrams. The treatment of 
non-perturbative effects in heavy quark scattering was first carried out by Braaten et al. [26-28] in thermal perturbation theory. 
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denoted as hard thermal loop (HTL) approach and later by Peshier et al. [29-31] and Gossiaux et al. [32-38]. In Ref. [39] 
we considered all the effects of the non-perturbative nature of the strongly interaction quark-gluon plasma (sQGP) constituents, 
i.e. the large coupling, the multiple scattering etc., where we refrain from a fixed-order thermal loop calculation relying on 
perturbative self-energies (calculated in the limit of infinite temperature) to fix the in-medium masses of the quarks and gluons 
and pursue instead a more phenomenological approach. The multiple strong interactions of quarks and gluons in the sQGP 
are encoded in their effective propagators with broad spectral functions. The effective propagators, which can be interpreted as 
resummed propagators in a hot and dense QCD environment, have been extracted from lattice data in the scope of the DQPM 
[22, 40,41]. 

In this work we extend our study of the scattering processes qQ qQ and gQ gQ in Refs. [23, 39] for a partonic medium 
at finite temperature T and in particular at finite chemical potential jlq. The gluons (g) and light (q) and heavy (Q) quark masses 
as well as the fundamental ingredients (infrared regulator TR’ and running coupling ’a^’) involved in the scattering amplitudes 
are determined in the framework of the DQPM. The dependencies of these quantities on T and jlq will be discussed below. 


A. Dynamical Quasi-Particles at finite temperature and quark chemical potential fj^q 

The DQPM describes QCD properties in terms of single-particle Green’s functions (in the sense of a two-particle irreducible 
(2PI) approach) and leads to the notion of the constituents of the sQGP being strongly interacting massive effective quasi-particles 
with broad spectral functions (due to the high interaction rates). The strategy for the determination of parton masses and widths 
within the DQPM approach is to fit the analytical expression of the dynamical quasi-particle entropy density to the IQCD 
entropy density ” which allows to fix the few parameters present in by lattice data in equilibrium [20-22, 41]. 

The variation of parton masses as a function of the medium properties is described by the spectral functions which are (except 
of a factor) identical to the imaginary part of the retarded propagator. These are no longer 5— functions in the invariant mass 
squared (as in the case for bare masses) [40]. For the current analysis, we use the approximation of momentum-independent 
real and imaginary parts of the retarded self-energy, which are - for a given temperature T - proportional to the parton mass 
and width, respectively [22]. A partonic propagator A is expressed in the Lehmann representation in terms of the spectral 
function A(;?) = / ff An often used ansatz to model a non-zero width is obtained by replacing the free spectral function 

Aq{p) = 27t[5{(0 — — 5(m + p^Y] by a Lorentzian form [22, 39] which also emerges from Kadanoff-Baym theory in first 

order gradient expansion. 

The extension of the DQPM to finite quark chemical potential Pq is more delicate since a guidance by IQCD is presently 
very limited. In the simple quasiparticle model one may use the stationarity of the thermodynamic potential with respect to 
self-energies and (by employing Maxwell relations) derive a partial differential equation for the coupling g^{T^Pq) which may 
be solved with a suitable boundary condition for g^{T^Pq = 0) [42-45]. Once g^{T,Pq) is known one can evaluate the changes 
in the quasiparticle masses with respect to T and /i^, i.e. dM^jdpq and dMljdT (for x = g^q^q) and calculate the change in 
the ’bag pressure’ (cf. Refs. [42-46] for details). However, such a strategy cannot be taken over directly in the DQPM since 
additionally the quasiparticle widths Yx{T,Pq) have to be known in the (T^Pq) plane in this case. 

In hard-thermal-loop (HTL) approaches [47, 48] the damping of a light quark (or gluon) depends weakly on the quark chemical 
potential explicitly [49]. This, however, has to be considered with care since HTL approaches primarily address Landau damping 
and assume small couplings g^. Accordingly, these concepts should be applied at sufficiently high temperature, only. Present 
IQCD calculations suggest that the ratio of pressure to energy density, P/e, is approximately independent of Pq as a function of 
the energy density e [50]. The functional dependence of the quasiparticle width 7 on T and Pq thus has to be modeled in line 
with ’lattice phenomenology’. We proceed in this paper to some scaling hypothesis in order to extend the definition of DQPM 
masses and widths to a finite chemical potential [22]. Assuming three light flavors (q = u^d^s) and all chemical potentials to be 
equal (Pu = Pd = Ps = Pq) the gluon and light quark masses are taken for T and Pq as 


This functional form is inspired by HTL masses but slightly differs from those 
the effective temperature for Nf = Nc = 3 



(II.l) 

in order to incorporate an explicit scaling with 
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which implies no additional parameter in the DQPM. Thus Mg{T^ fiq) and Mq{T^ jiq) are given by 




g {T /Tci^lg)) ^ ^* 2 ^ 


(II.3) 


with the coupling constant g^{T'^/Tc{llq)) in (11.3) is considered here as depending on the medium temperature and for > Tg 
is given by: 


g\T\^,) = 


48;r^ 


{nNc-2Nf)ln (^A2( 


T'^ T 

Tcil^q) Tc{lXq) 


?) 


(II.4) 


Since the coupling (squared) in the DQPM is a function of T jTc a straight forward extension of the DQPM to finite jlq is to 
consider the coupling as a function of T""lTc{llq) with a /i^-dependent critical temperature, 


= ■■ 

with a 8.79 GeV“^. The expression of Tc{llq) in (II.5) is obtained by requiring a constant energy density £ for the system 
at r = Tc{llq) where £ at Tc{llq = 0) « 0.158 GeV is fixed by lattice QCD calculation at jlq = 0. The coefficient in front of the 
-dependent part at first sight appears arbitrary but can be compared to recent IQCD calculations (for imaginary quark chemical 
potentials) at finite (but small) {ds which gives [51] 


Tc{I^b) ^ 
Tc 
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with K = 0.013(2). Rewriting (II.5) in the form (II.6) and using fis ~ Sjlq we get Koqpm ~ 0.0122 which compares very well with 
the IQCD result. Consequently one has to expect an approximate scaling of the DQPM results if the partonic width is assumed 
to have the form, 
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This choice leads to an approximate independence of the potential energies per degree of freedom as a function of fdq. Never¬ 
theless, the conjecture (II. 7) should be explicitly controlled by future IQCD studies for Nf=3 at finite quark chemical potential. 
Unfortunately, this task is presently out of reach and one has to live with the uncertainty in (II. 7) which is assumed in the following 
investigations. 

We recall that within the scaling hypothesis (II.2)-(II.7) the results for the masses and widths at finite T stay about the same as 
a function of T""/Tc{llq) when dividing by the temperature T [22]. 

For a finite quark chemical potential jlq the energy density £ in the DQPM is seen to scale well with {T/Tc{llq = ff)Y as a 
function of temperature for T""lTc{llq) > 3, however, increases slightly with jlq close to the phase boundary where the scaling 
is violated on the level of 20%. This violation in the scaling is essentially due to an increase of the pressure P [22]. Note that 
a quark chemical potential of 0.21 GeV corresponds to a baryon chemical potential of flp ^3flq =0.63 GeV which is already 
substantial and the validity of (II.5) becomes questionable. Since the pressure P is obtained from an integration of the entropy 
density 5' over temperature T the increase in P with jlq can directly be traced back to a corresponding increase in entropy density. 

Next we discuss the influence of a finite chemical potential jlq on the running coupling a^, the DQPM masses, widths and 
parton spectral functions. The running coupling /{An) is presented in Fig. l-(a) as a function of the temperature T 

for Jlq = f) and in Fig. l-(b) for finite T and jlq. One sees that is larger than 1 near Tc{llq) and non-perturbative effects 
are most pronounced at these temperatures, with Tc{liq = 0) = 0.158 GeV and Tc{liq = 0.2GeV) 0.127 GeV. Note that close 

to Tc{llq = 0) the full coupling calculated on the lattice increases with decreasing temperature much faster than in the pQCD 
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EIG. 1. (Color online) (left) The DQPM running coupling as{T,Pq) = g^ (T'*'/Tc{pq))/{47l) as a function ofT in the IQCDfor Nf = 0 (red 
spheres) [52] and in the DQPM at zero chemical potential for Nf = 0 (black line) and Nf = 3 (dashed brown line). The result for the DQPM 
at finite quark chemical potential (Pq = 0.2 GeV) for Nf = 3 is given by the orange thin line, (right) 3D plot ofas(T,Pq). 




EIG. 2. (Color online) The DQPM pole masses and widths for the gluons (Mg, Yg) (left) and light quarks (Mq, jq) (right) given by (II.l) and 
(IL7) as a function of temperature for different values of the quark chemical potential Pq. 


regime (at high T). A finite Pq leads to a smaller value of the coupling constant as compared to the Pq = 0, except close to the 
corresponding critical temperature Tc{Pq). Indeed, the value of as{Tc{Pq = Q)) =2.84 while as(Tc(Pq = 0.2 GeV)) =4.7. For 
temperatures larger than Tc{pq = 0) one finds systematically as{T^Pq = 0) > as{T,Pq), as shown in Fig. l-(b). 

The DQPM masses and widths for the gluon (g) and light quarks (q), given by (II.l) and (II.7), are presented in Fig. 2-(a) 
and (b), respectively. Since the width is found to be much smaller as the pole mass, the excitations can well be considered as 
dynamical (off-shell) quasiparticles. For larger T, after a shallow minimum at T ^ 1.2Tc{Pq), the width / increases slowly with 
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T and even for large T is to a good accuracy proportional to the temperature (and also to the mass). 

In Figs. 2-(a) and (b), the masses and widths are plotted as a function of the temperature. It is seen that the gluon and light 
quark masses and widths decrease at fixed T with increasing quark chemical potential jlq. The finite flq has a larger effect on 
the masses as compared to the widths. For completeness we note: Tc{liq = 0) = 0.158 GeV, Tc{liq = 0.2GeV) = 0.127 GeV and 
Tcildq = O.SGeV) = 0.072 GeV. 

Using the pole masses and widths (II. 1) and (II. 7) and the DQPM running coupling (cf. Fig.l), the Breit-Wigner spectral func¬ 
tions for the different QGP species are completely determined, i.e. the imaginary parts of the retarded propagators. Accordingly, 
also the full retarded propagators for the effective partonic degrees of freedom are known in first order gradient expansion and 
we can proceed with the calculation of some amplitudes keeping in mind the uncertainties at high fdq. 


III. qQ AND gQ ELASTIC SCATTERING AT FINITE T AND giq 

The matrix elements for the qQ qQ and gQ gQ channels have been calculated for the case of massless partons in the 
vacuum in Refs. [24, 25]. In this section we study the qQ and gQ elastic scattering in the QGP medium at finite temperature T 
and chemical potential jlq considering the case of on- and off-shell gluons and light and heavy quarks. The partons are dressed by 
effective masses in the on-shell case and are dressed by the DQPM spectral functions with a finite width in the off-shell case. We 
refer to the on-shell study as DpQCD approach (Dressed pQCD) and to the off-shell case by lEHTL approach (Infrared Enhanced 
Hard Thermal Loop). For the DpQCD approach the qQ and gQ elastic cross section at finite T and jlq is determined by using (i) 
the running coupling as{T^iiq) (Fig. 1), and (ii) the DQPM pole masses for the incoming and outgoing quarks and gluons. The 
DQPM gluon pole mass serves also as an infrared regulator in the gluon propagator. 

Considering in- and out- dynamical quasi-particles (DQP), the corresponding quasi-elastic lEHTL cross section for the 

process (1)"^^^^ + (2)"^^^^ ^ (3)"^^^^ -h (4)"^^"^^ is deduced by the convolution of the modified pQCD cross section CT, where complex 
propagators are considered in the transition matrix elements, with the spectral functions, i.e. 

^IEHTL(^) = jdmW pf^{m^^^), (ULl) 

where p^.^(m^‘^) is the Breit-Wigner spectral function of the particle i, normalized as J dm^'^p ^^= 1. The on- and 

off-shell cross sections for qQ, gQ scattering in the sQGP at finite T and jlq are obtained using the DQPM parametrizations for 
the quark (gluon) self-energies, spectral functions and interaction strength at finite T and flq. Here, in principle, a two-particle 
correlator should appear, but since we work in a 2PI motivated scheme the partons in the sQGP can be characterized by (dressed) 
single-particle propagators (cf. Ref. [23]). 

In the context of the hot and dense QGP, the elementary Feynman diagrams for the qQ and gQ elastic scattering at order 0{as) 
are illustrated in Fig. 3. 



FIG. 3. (Color online) Feynman diagrams for the qQ qQ and gQ —)■ gQ scattering process. Latin (Greek) subscripts denote colour (spin) 
indices, ki, resp pt (kf, resp pf) denote the initial (final) 4-momentum of the light quark or the gluon, resp the heavy quark. The invariant 
energy squared is given by s = (pt -h kt)^, t = (pt — pf)'^, u = (pt — kf)'^. 


A. On- and off-shell qQ elastic scattering 

The process qQ qQ is calculated here to lowest order in the perturbation expansion using the extended Feynman rules 
for massless quarks in Politzer’s review [53] for the case of finite masses and widths. The color sums are evaluated using the 
techniques discussed in Ref. [53]; the spin sums will be discussed below. Contrary to the case of massless gluons where the 
'‘Transverse gauge ” is used, the "Lorentz covariance ” is used for the case of massive gluons here since a finite mass in the gluon 
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propagator allows to fix the O’th components of the gluon fields (a = 1, • • • , 8) by the spatial degrees of freedom A^ (^ = 1,2,3). 
Furthermore, the divergence encountered in the Uchannel (Ref. [24, 25, 54]) - when calculating the total cross sections and 
G^^ - is cured self-consistently in the DpQCD and lEHTL models since the infrared regulator is given by the finite DQPM gluon 
mass (and width) in the DpQCD (lEHTL) model. For on-shell qQ elastic scattering, the Uchannel invariant squared amplitude - 
averaged over the initial spin and color degrees of freedom and summed over the final state spin and color - Mt is given by 


IIMP 


9(t-m^Y 


+ {u-M^-mlf+ 2{hfQ + ml)t , 


(IIL2) 


where ruq (Mq) is the light quark (heavy quark) mass and rUg is the DQPM exchanged gluon mass. 

In the off-shell picture we take into account not only the finite masses of the partons, but also their spectral functions, i.e. their 
finite widths. Since the light quark and heavy quark masses change before and after the scattering (’quasi-elastic’ process) we 
introduce the mass for the initial q and for the final q, and allow for different masses of the heavy quark, Mq for the initial 

Q and Mq for the final Q. The squared amplitude - averaged over the initial spin and color degrees of freedom and summed over 
the final state spin and color - gives: 






9[(r-m|)2 + 47|^2j _ 


[pfPj +pfp} + g^''!^) ‘^(kf,^ki^v+ki,^kf^y+g^y^-'j 


(IIL3) 


where we have incorporated the DQPM propagators (i.e. = t — m^± ^ijgqo, with rUg, Jg is, respectively, the effective gluon 

mass and total width at temperature T and quark chemical potential jlq and — k? is the gluon energy in the 

^-channel). Thus the divergence in the gluon propagator in the t channel is regularized. 

The relative contribution of the off-shell partons to the pQCD cross section is expected to change due to different kinematical 
thresholds and to the changes in the matrix element- corresponding to the diagram in Fig. 3. The off-shell kinematical limits for 
the momentum transfer squared t and the expressions of the Mandelstam variables in the case of off-shell heavy quark scattering 
are given in Ref. [39] . 


B. On- and off-shell gQ elastic scattering 

The invariant amplitudes for the three graphs (shown in Fig. 3) for the case of massive heavy quarks and massless gluons is 
given according to Combridge [25] and revisited in [39]. In Ref. [39] we have already studied the scattering of massive on- and 
off- shell heavy quarks on massive on- and off-shell gluon in a finite temperature medium at jlq = 0. In this paper we extend this 
study to a medium at finite temperature and chemical potential. Therefore, the lowest-order amplitude for the process gQ gQ, 
obtained from the Feynman rules of the gauge theory by the sum of the amplitudes of the three graphs (cf. appendix A.4 of [39]), 
where in order to obtain the correct result for the squared matrix element 

< \M\^ >= I L TapT*,p,s^e*‘^'efsf, (IIL4) 

^ spins 

we have to use appropriate projection operators for the transverse polarisation states 


^i,oc^i,a' — Saa' 

pol,i 


^i,oc ^i^a' 


E =gpli'- 

polj 




(III.5) 


We recall that for vector fields with nonzero Lagrangian mass there is no gauge freedom anymore. The massive vector field 
A^ only has to fulfil the condition 5^A^ = 0. Therefore, the sum over the initial and final gluon polarizations is fixed by the 
expressions in (III. 5). 

For the case of finite masses and widths of the scattering quarks and gluons, the gQ elastic scattering amplitude has been given 
in Ref. [39], where we have to take into account the spectral functions for the heavy quark and gluon masses at finite temperature 
and chemical potential, the coupling (cf. Fig.l) at finite T and jlq and the quark and gluon propagators for the case of massive 
vector gluons with finite lifetime {q^ nig) and for the case of massive fermions with finite life time Sfip^niq): 
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GY\q) = -i 




H 


/-m2 + /27g(pi)-p^)' 


^f{p) — 


f + MQ 


u-Ml + i2YQ{pi^-k{y 


^f{p) — 


f + MQ 


s-M^ + ilygipi^+pl)' 


(III.6) 


where rUg, Jg (Mq, Jq) are the mass and width of the gluon or the heavy quark at finite temperature and chemical potential (cf. 
Sec.II). We note that the heat bath breaks the Lorentz covariance. Accordingly the energies and denote quantities in the 
rest frame of the heat bath. 


C. Results for elastic scatterings 

We consider first the scattering of a (high momentum) heavy quark with a light quark or a gluon in a QGP at temperature T 
= 0.2 GeV with invariant energy = 4 GeV for different quark chemical potential Pq = 0,0.2,0.3 GeV. Fig. 4-(a) presents the 
off-shell differential cross section da/d cos 6 (black lines) in comparison to the on-shell cross section (orange lines) of the uc 
elastic scattering. The importance of finite width corrections in the uc scattering processes is illustrated by comparing the two 
differential cross sections. For the energy of = 4 GeV one observes a deviation of the off-shell results compared to the on- 
shell ones only for large scattering angles. However, according to the small differences between the differential on- and off-shell 
cross sections one can conclude that the total on-shell cross sections do not change on a relevant scale when introducing off-shell 
masses. This is particularly important since the width of the heavy quark has been taken as an upper limit (cf. Section II). Figure 
4- (a) shows also the influence of a finite chemical potential on the heavy quark scattering. The increase of Pq leads to a decrease 
of da/d cos 0 and consequently to a decrease of the total cross section as illustrated in Fig.4-(b). Despite the lower value of the 
IR regulator at higher values of Pq (decrease of the gluon mass at finite Pq) the smaller values of the running coupling at finite Pq 
explain the decrease of da/d cos 6 at finite and large values of Pq. 



FIG. 4. (Color online) Differential (a) and total (b) elastic cross section for uc uc elastic scattering for off-shell (black lines) and on-shell 
partons (orange lines) at three different values of the quark chemical potential Pq (see legend) at T = 0.2 GeV We consider the DQPM pole 
masses for the on-shell partons and the DQPM spectral functions for the off-shell degrees of freedom. 

The total elastic cross section of a c-quark, which traverses a plasma at temperatures T = 0.2 GeV, as calculated in the DpQCD 
and lEHTL approaches, is shown in Fig. 4-(b) as a function of y(y for different values of Pq. Apart for energies close to 
the threshold the cross sections show a rather smooth dependence on the invariant energy yT, however, differ substantially in 
magnitude with temperature and chemical potential. Fig. 4-(b) demonstrates also that, independently on Pq, the off-shell mass 
distributions only have a sizeable impact at the threshold given by the pole masses for uc scattering. This is due to the moderate 
parton widths considered in the DQPM model. At energies below the on-shell threshold the off-shell cross section increases with 
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^/s because more and more masses can contribute. Whereas the on-shell cross section diverges at the threshold the off-shell cross 
section shows a maximum at the on-shell threshold and decreases then due to the decrease of the on-shell cross section. 

The conclusions drawn before for the study of uc elastic scattering are valid also for gc scattering, however, with cross sections 
for gQ elastic scattering that are larger than the cross sections for uc scattering by roughly a factor of 9/4 which is ratio of 
the different color Casimir operators (squared). This is also related to the fact that the scattering of heavy quarks with gluons 
proceeds Vmt,s and u channels, whereas one has only the t channel for uc elastic scattering. 

Figs.5-(a) and (b) show explicitly the temperature and ^/s dependences of the uc and gc elastic cross sections at jlq = 0, as 
described in the DpQCD approach. We deduce that an increasing medium temperature T leads to an increase of the thermal 
gluon mass (infrared regulator) and hence to a decrease of the DpQCD uc and gc elastic cross sections. The effective gluon mass 
is roughly proportional to T for temperatures above 0.2 GeV. The large enhancement of the total cross section for temperatures 
close to Tc{llq) can be traced back to the infrared enhanced coupling in the DpQCD/IEHTL models. 



EIG. 5. (Color online) Elastic cross section of uc ^ uc (a) and gc ^ gc (b) scattering as a function of the temperature T and the invariant 
energy above threshold y/s — where is the threshold energy, for on-shell partons as described by the DpQCD approach at Pq = 0. 

In order to quantify the effect of a finite Pq on the heavy quark scattering, we show in Figs. 6-(a) and (b) the uc and gc thermal 
transition rate CO (defined by Eq. (3.8) in Ref. [23]) as a function of T and Pq. We provide the results only for DpQCD because 
there is just a small difference between the on-shell DpQCD and the off-shell lEHLT approaches (cf. Fig.4). 

There are two different profiles in the {T,Pq) dependencies of the thermal transition rate CO. For temperatures larger than 
Tc{Pq = 0) =0.158 GeV, a small effect of a finite Pq, leading to a decrease of co, is noticed. For temperatures smaller than 
Tc{Pq = 0), an increase of CO appears when Pq increases. The last effect is due to the increase of the coupling at these 
temperatures. More precisely, one has as{T < Tc{Pq = 0),/i^) > CXs{T = Tc{Pq = 0),Pq). Besides the transition amplitudes, one 
should also consider the effect of the statistical weights (Fermi-Dirac distribution function) at finite (T,Pq) in the evaluation 
of CO; /™ increases for large values of Pq leading to an extra contribution to the increase of the thermal transition rate CO for 
T < Tc{Pq = 0), but is not enough to counterbalance the decrease of the total cross section for T > Tc{Pq = 0). 

The thermal transition rate can be parametrized by a power law in T for each value of Pq, i.e. ^ T~^ for T > Tc{pq = 0). In 
fact, one can find that p^c{Pq=o)<T<i.2Tc{pg) ^ ^ ijT>i.2Tc{pq) ^ 2 . These different power laws in T will have a sizeable effect on 
the transport coefficients to be evaluated in the following sections. On the other hand, having almost the same power laws at finite 
Pq as compared io Pq = Q (for temperatures larger than Tc{pq = 0)) leads to some scaling effects in the transport coefficients at 
these temperatures. 

Additionally, it is worth to evaluate the uc and gc elastic scattering cross section as a function of the energy density £ available 
in the heavy-ion collisions since out-off equilibrium a temperature is ill defined while £ can be well calculated in the local rest 
frame (e.g. in PHSD). Here the energy density £ for a given temperature T and quark chemical potential Pq is obtained by using 
the inverted DQPM equation of state which gives the temperature as a function of the energy density £ for a given quark chemical 
potential. The temperature as a function of £ for = 0,0.1,0.2 GeV is shown in Fig. 7-(a). Note that the DQPM model describes 


9 










H. BERREHRAH et al 



EIG. 6. (Color online) Thermal transition rate CO of uc ^ uc (a) and gc ^ gc (b) as a function of the temperature T and quark chemical 
potential Pqfor on-shell partons as described by the DpQCD approach. 


the QCD energy density at temperatures even as low as T ^Tc{liq). 



e [GeV/fm^] € [GeV/fm^] 


EIG. 7. (Color online) Temperature T as a function of the energy density £ from the DQPMfor different values of the quark chemical potential 
Pq (a); Thermal transtion rate for uc uc scattering for off-shell (lEHTL) and on-shell (DpQCD) partons as a function of the energy density 
8(b). 


Transport theoretical simulations have shown that the local energy densities achieved in the course of heavy-ion collisions at 
FAIR energies reach at most 3 GeV fm“^, at SPS and RHIC energies up to 30 GeV fm“^ and up to 300 GeV fm“^ at LHC. 
Therefore one observes that the uc elastic cross section, displayed in Fig. 7-(b) at the energy densities of interest is G [0.08-8] mb 
from LHC to FAIR energies, following the DpQCD/IEHTL approaches, with typical values of 5 mb at the phase transition line. 


10 













COLLISIONAL PROCESSES OE HEAVY QUARKS ... 


IV. HEAVY QUARK INTERACTION RATES IN A MEDIUM AT FINITE T AND 


Using the elastic cross section for q{q)Q and gQ collisions, for on- and off-shell partons - as calculated in Sec.III - we evaluate 
the interaction rate of a heavy quark with momentum p and energy E propagating through a QGP in thermal and chemical 
equilibrium at a given temperature T and quark chemical potential jlq. The quarks/antiquarks of the plasma are described by a 
Fermi-Dirac distribution fq^q{q) = {EqT^q)/T^^ whereas the gluons follow a Bose-Einstein distribution fg{q) = ^ • 

For on-shell particles (DpQCD model) and in the reference system in which the heavy quark has the velocity P = pjE the 
(on-shell) interaction rate R'^^{p) 


for 2 —2 collisions is given by [23] 


dt 




r^'{s)fr{q) 


sEn 


dq^ 


(IV.l) 


where ^ denotes the sum over the light quarks/antiquarks and gluons of the medium. In Eq. (IV.l) fr{q) is the invariant 
distribution of the plasma constituents in the rest frame of the heavy quark, given for the quark/antiquark by: 


J dO. fr{q) =2n J dcosGr 


U qCOsQr^l^q)/T _j_ 


(IV.2) 


with u = (iP^u) = ^{E^ —p) being the fluid 4-velocity measured in the heavy-quark rest frame, while Or is the angle between 
q and u. m^{s) in (IV.l) is related to the transition amplitude |A^2,2p of the collision q{q^g)Q q{q^g)Q by 




1 1 


2pLi^)J-4pL SQgpiJkJ 




(IV.3) 


with pcm = {q Mq) / ^/s denoting the momentum of the scattering partners in the c.m. frame and gQ (gp) the degeneracy factor 
of the heavy quark (parton). 

For off-shell partons (lEHTL model) the elastic interaction rate is obtained by replacing 


/ 


(fip 1 
{2nf W ' J {2n) 


-A 


/Ap(/>)e(M) 


(IV.4) 


with p{p) denoting the spectral function which can be specific for each particle species. Eor the Breit-Wigner-m form of the 
spectral function, the explicit extension of Eq. (IV.l) for the off-shell case is given by 


R°^{p,T,Pq) = E ^iep,q,p',q' J {nii) J 


q^mf{s)friq) 


s E„ 


dq, 


(IV.5) 


where 


mf{s) = 


2p^p 


\ rtm 
Juir. 


Sp§Q ij 


(IV.6) 


with Y,\M 2 ^ 2 iP^^dJ\p',q'\kJ)\^{s,t) being the off-shell transition amplitude defined in the lEHTL approach (cf. Ref.[39]); 
p^ {pf) is the initial (final) heavy-quark momentum and the integration boundaries tmin and tmax are given by: 


= -| (^1 - (ft +ft + ft +ft) + (ft -ft)(ft - ft) ± ^(l-ft-ft)2-4ftft^(l-ft-ft)2-4ftftj , 

with: ft = {m\f Is, ^ (IV.7) 
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The total interaction rate (IV. 5) of the off-shell approach (lEHTL) in the plasma rest system is compared to that of the on-shell 
calculations (DpQCD) (IV. 1) in Fig. 8-(a) as a function of the momentum of the heavy quark p for different values of the quark 
chemical potential jlqdiiT = 0.2 GeV. We assume here a Breit-Wigner spectral function and a Boltzmann-Jutner distribution for 
both, the light quarks/antiquarks and the gluons. Our results are rather independent of the choice of the spectral function or by 
replacing the Boltzmann-Jutner distribution by a Fermi/Bose distribution. In Fig. 8-(a) the black lines refer to lEHTF results and 
the orange lines to DpQCD results. We see from Fig. 8 -(a) that the finite width of the spectral function (with the DQPM width) 
decreases the interaction rate of heavy quarks with the medium on the order of 20%. This modification is rather independent of 
the heavy-quark momentum, temperature and quark chemical potential of the plasma. The difference between the DpQCD and 
lEHTF rates is related on one side to the propagator, which contains an additional imaginary part proportional to the gluon width 
in the lEHTF model, and on the other side to the energy asymmetric contribution of the Breit-Wigner spectral function in lEHTF. 
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FIG. 8. (Color online) The total elastic interaction rate R of c-quarks in the plasma rest frame in the lEHTL and DpQCD models as a function 
of the heavy-quark momentum pfor three different values of the quark chemical potential Pq = 0,0.2,0.3 GeV atT = 0.3 GeV (a), (b) Contour 
plot of the total elastic interaction rate R within DpQCD as a function ofT and Pq. 


For fixed temperature T = 0.3 GeV the variation of the quark chemical potential Pq from Pq = 0 leads to a decrease of the 
rate (cf. Fig. 8 (a)). However, e.g. for the temperature T = 0.2 GeV an increasing Pq leads to either an increase of the rate (for 
Pq = 0.2 GeV) or to the decrease of R (for Pq = 0.3 GeV). Therefore, the dependence of R on both the medium temperature and 
quark chemical potential is not trivial. The dependence of the rates on the medium temperature T and quark chemical potential 
Pq for a heavy-quark with a momentum p = 5 GeV is illustrated in Fig. 8-(b) in DpQCD. As expected, the rate increases for 
higher temperatures dX Pq = 0 because the number of plasma particles becomes larger T^). Therefore, the increase of the 
interaction rate with temperature keeping Pq small is seen. For a fixed temperature, the variation of Pq leads to different profiles 
in the rates. A decrease of the total rate for high temperatures when Pq increases and an opposite trend for small temperatures 
(T <Tc{Pq = 0))\^ observed. The highest values of the rates are reached for small Pq and large temperatures T 

Due to the different abundances of particle species in a medium at finite chemical potential, it is interesting to study the variation 
of the heavy-quark interaction rates with the quarks/antiquarks and gluons independently. Figs. 9-(a), (b) and (c) illustrate the 
dependence of the heavy-quark collisional rates with quarks, antiquarks and gluons of a medium at finite temperature T and quark 
chemical potential Pq for an intermediate heavy-quark momentum (p = 5 GeV). 

For the case of gluons and antiquarks, the interaction rates are increasing with higher temperature for all Pq. The charm quark 
interaction rate with light quarks (Rue) depends on (T,Pq) in a similar fashion as described in Fig. 8-(b) for large Pq and small 
temperatures (T <Tc(Pq = 0)). Indeed, for larger values of Pq and small temperatures. Rue is much larger than Rue and Rgc, so 
that the total interaction rates are dominated by Rue^ On the other hand the R profile is dominated by Rgc for small Pq and large 
temperatures. This is easy to interpret: At large Pq the number of light quarks is large compared to the number of antiquarks, i.e. 
the Fermi-Dirac distribution contributes unequally to R for u and u. On the other hand the gluon number decreases with larger Pq 
since it is correlated with the subdominant light antiquarks, via the T and Pq dependencies of the masses. 
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EIG. 9. (Color online) The total elastic interaction rate R of c-quarks in the plasma rest frame as a function of the temperature T and quark 
chemical potential Pq due to the scattering with light quarks (a), light antiquarks (b) and gluons (c). The on-shell heavy-quark momentum in 
all cases is p = 5 GeV 


V. HEAVY QUARK ENERGY LOSS IN A MEDIUM AT FINITE TEMPERATURE T AND CHEMICAL POTENTIAL Pq 

The collisional energy loss dE/dt has been formulated by Bjorken and been explicitly calculated in Ref. [23] for Pq = 0. We 
recall here the expression of dEjdt for on- and off-shell partons. In the framework of the DpQCD model, dEjdt is given in the 
plasma rest frame by: 


dE^^ 

dt 


{P,T,flq) 


Mq P 
4(27r)3 Jo s^Eq 


-qfoiq) + ^ {Mq + Eq)fi (q) 


dq. 


Eq.(V.l) is easily extended to the off-shell case with Breit-Wigner spectral functions using (IV.4). One finds [23] 


(V.l) 


dE°^. 
dt 


J Pi if^i) 


(s) 

S^En 


-qh{q) + ^{MQ+Eq)p{q) 


dq, (V.2) 


where m“(s), m®" (s) and /i {q) in (V.l) and (V.2) are given by 

^Pcm SpSQ ij P I 


„„ 1 ptmax 1 ___ ___ ^ 

/ ——Y.T.\^f2{s,f,iJ\k,l)\ 

^P'P^ Jtmin gpSQijtl 

1 /'j a r f Eq-uq cos Or TPq\ no 

fn{q) = 27 dcosdr f ( - - - - - - 1 cos" dr . 


1 - 


t - {M'qP + {M^qP - 2E‘Er 
2p^pf 


dt, 


(V.3) 


Furthermore, p^ (Mq) is the initial heavy momentum (mass) and p^ (Mq) is the final heavy momentum (mass). The heavy- 
quark energy loss (eq.V.l) and (eq.V.2) is illustrated in Fig. 10 as a function of the heavy quark momentum for the DpQCD and 
lEHTL models at T = 0.2 GeV for different quark chemical potentials up to Pq = 0.3 GeV. As for the rates, the off-shell spectral 
function decreases the energy loss as compared to the on-shell case as described by the DpQCD model. The difference is uniform 
as a function of the heavy quark momentum p and the increase of the quark chemical potential. For low heavy-quark momenta 
and higher medium temperatures, the heavy quark looses less energy by elastic collisions. At low momentum heavy quarks start 
to gain energy to approach thermal equilibrium and to arrive at the average energy in the heat bath (dE/dx is negative). For a 
medium at the temperature T = 0.2 GeV, the increase of Pq leads to a decrease of the heavy-quark energy loss. Nevertheless, 
varying both the medium temperature and Pq influences the heavy-quark energy loss in different ways. 
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EIG. 10. (Color online) The c-quark energy loss, dE/dx, 
in the plasma rest frame from the DpQCD and lEHTL 
models as a function of the heavy-quark momentum p at 
T = 0.2 GeV for different values of the quark chemical 
potential Pq. 


The temperature and quark chemical potential dependencies of the collisional energy loss for low (p = I GeV/c) and interme¬ 
diate heavy-quark momentum (p = 5 GeV/c) is displayed in Fig. 11 -(a) and (b) for the on-shell DpQCD model, respectively. 
Figs. 11- (a) and (b) show that the energy loss is increasing as a function of p. For intermediate momenta p the conclusions 
drawn for the rate are valid also for the collisional energy loss; especially Fig. 1 l-(b) shows that with increasing temperature for 
small values of Pq - which corresponds to the increase of the parton masses in the DpQCD- the energy loss becomes larger. This 
is more pronounced for large heavy-quark momenta. On the contrary, for p = 1 GeV (cf. Fig 11 -a), where we observe an energy 
gain (a negative energy loss), the energy gain increases with the distance of the energy of the heavy quark from its equilibrium 
value. Nevertheless, as expected for a cross-over transition (at finite but small values of Pq) we observe a very smooth dependence 
of the energy loss on both variables, T and Pq. The dependencies of dE/dx on T and Pq are studied in [49] using an extension of 
finite temperature HTL calculation to include finite Pq in the distribution functions of the quark/anti-quarks. The effect of finite 
Pq on dE/dx seen in [49] is much smaller than in our DpQCD and lEHTL models, shown in Fig. 10. We note that the finite 
Pq effects in both the parton masses, coupling constant and distribution function are included in our model, which give to our 
calculations more inclusive effects than those of [49]. 

The energy loss of a heavy quark due to its scattering with the QGP partons varies with the abundance of the particle species. 
Whereas the scattering of heavy quarks on the gluons is slightly Pq dependent, the scattering on light quarks/antiquarks induces 
an energy loss which is highly (T, Pq) dependent. Figs. 12-(a), (b) and (c) illustrate this dependence of the heavy-quark 
collisional energy loss with quarks, anti-quarks and gluons of a finite temperature T and quark chemical potential Pq medium 
for an intermediate heavy-quark momentum (p = 5 GeV). As seen for the interaction rate, the profiles observed for the q{q) — Q 
collisional energy loss evolve oppositely due to the Fermi-Dirac distribution. As pointed out before, the quark (antiquark) density 
changes substantially with Pq. 


VI. HEAVY QUARK MOMENTUM LOSS IN A MEDIUM AT FINITE TEMPERATURE AND CHEMICAL POTENTIAL 

We continue our study with the momentum loss of heavy quarks in a partonic medium at finite temperature and chemical 
potential. The drag (A) and diffusion (B) coefficients are evaluated according to [32, 55, 56] by a Kramers-Moyal power expansion 
of the collision integral kernel of the Boltzmann equation. Note that the diffusion tensor B admits a transverse-longitudinal 
decomposition (perpendicular and along the direction of the heavy quark in the fiuid rest frame) and contains two independent 
coefficient 5|| and 5^. 


A. Longitudinal momentum loss 


The drag coefficient A describes the time evolution of the average of the longitudinal component of the momentum transfer 
{p ~ p'Y of the heavy quark. For on-shell partons it has been defined in the plasma rest frame by Svetitsky [55, 56] and evaluated 
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EIG. 11. (Color online) The c-quark energy loss, dE/dx, in the plasma rest frame from the DpQCD approach as a function of the temperature 
T and quark chemical potential Pqfor the heavy-quark momentum p = 1 GeV (a) and p = 5 GeV (b). 



EIG. 12. (Color online) The energy loss dE/dx ofc-quarks in the plasma rest frame as a function of the temperature T and the quark chemical 
potential Pq due to the scattering with light quarks (a), light antiquarks (b) and gluons (c). The on-shell heavy quark momentum in this case is 
p = 5 GeV 


in [23] by 


A^yp,T,^g) = 


d<p> Mq 


■[ 




dt A{2 kY Jo s^E^ 
Eq.(VI.l) is easily extended for the off-shell case. One finds [23], 


{MQ+Eq)f,{q)-Uh{q) 


dq. 


(VI.l) 


ni^p,q,p’,q' J mpiitidmi pf^irtii) / 


q‘^m^{s) 

s'^Ea 


{MQ-^Eq)fi{q) - ^qfo{q) 


dq, (VI.2) 
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where and /o,i {q) are given in (V.3). 


B. Transverse momentum loss 


During its propagation through the plasma a heavy quark receives random kicks (momentum transfers) from the constituents 
of the medium. The average value of the transverse momentum remains zero but the variance < p\> increases with the number 
of collisions. The increase of the variance <Pj> per unit length is called the transport coefficient q and defined as 




^ ^ P± ^ ^ P± ^single coll 


dx 


i 


(VI.3) 


Using the relation between the transverse second moment Bj and q given by ^ ^ we can evaluate q for massive on-shell 

partons (as in [23]) by the expression: 


r{pj.^q) = 


K 


L 


+ CX, ^5 


2{2k)^P Jo 


(A+A) + (A - A) 


dq, 


(VI.4) 


where m°"(s) and fo,i,2{q) are given in (V.3) and m‘^{s) by 




m^{s)= / dcosO 


l-cos 0 \^ 1 


SgSQ ij kl 


^^PcmJ-4pL SgSQ^a 


The off-shell transport coefficient q^^ is deduced by using Breit-Wigner spectral functions for the off-shell partons and extending 
the on-shell < 7 ®^ as 


T, Hq) = 3 ^iep,q,p',q'J —rtiidniipf^ {mi)J 


2(271) 


s^E„ 




(/0 + /2)- 


(Eq -\-mpY 


^2 W(/o-/2) 


dq (VI. 6 ) 


with: 


1 Ptmax 1 ___ ___ ^ 

= ——EEl-^ I 

^PP^Jt^in SggQ^U 


mf{s) = 


1 - 


t - (My 2 + (My 2 _ 2E‘Ef 

2p^pf 


dt. 


(VI.7) 


C. Numerical results 

We discuss now the drag coefficient and the transport coefficient q for the models DpQCD and lEHTL from Ref. [39] at 
finite T and Pq. In the DpQCD approach, the light and heavy quark and gluon masses are given by the DQPM pole masses. 
Fig. 13-(a) shows the drag coefficient (VI. 1) and (VI. 2) of heavy quarks as a function of the heavy quark momentum p. The 
temperature of the heath bath is chosen as T = 0.2 GeV. Figure 13-(b) illustrates the influence of a finite parton width on the heavy 
quark transport coefficient q, where q is displayed for on-shell partons (DpQCD) and off-shell partons (lEHTF) as a function of 
the heavy quark momentum and for a temperature of T = 0.2 GeV and different values of Pq. For both A and q the off-shell 
partons case is slightly lower by about 20 % independent of momentum, temperature or quark chemical potential than for the 
corresponding on-shell case. 

The temperature and quark chemical potential dependences of the drag coefficient and q for intermediate heavy quark momen¬ 
tum (p = 5 GeV/c) is displayed in Fig. 14-(a) and (b) for the on-shell DpQCD model, respectively. Fig. 14- (a) and (b) shows 
that both the longitudinal and transverse momentum is increasing for high temperatures and low Pq. On the contrary, for low 
temperatures and large Pq the longitudinal momentum loss is still considerable whereas the transverse momentum transfer is low. 
We conclude that longitudinal momentum transfers are important not only in a hot medium but also in a dense medium whereas 
the dense medium leads to less transverse fluctuations in the heavy quark propagation. The relatively large drag at low tempera¬ 
tures in DpQCD/IEHTF is due to the strong increase of the running coupling as{T^Pq) (infrared enhancement) for temperatures 
close X.oTc{Pq). 
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heavy quark momentum p for T = 2 GeV. 
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EIG. 14. (Color online) c-quark drag (a) and q (b) in the plasma rest frame from the DpQCD approach as a function of the temperature T and 
quark chemical potential Pqfor the heavy-quark momentum p = 5 GeV. 


Regarding the variations in T and jlq of the longitudinal and transverse momentum losses, one may roughly presage the ten¬ 
dency of the nuclear modification factor Raa and the elliptic fiow of charm particles in a hot and dense medium. We expect that the 
Raa of charmed mesons will still be small in a hot and dense medium as in case of finite temperature and zero chemical potential. 
On the contrary, the elliptic flow V 2 is expected to be much smaller in a hot and dense medium compared to our knowledge from 
RHIC and LHC measurements for jiq ^ 0, since the dense medium damps the transverse momentum fiuctuations. Moreover, one 
could notice that these finite jlq calculations might have interesting consequences on the c-quark angular correlations at FAIR, 
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which should be more peaked than those seen at RHIC or LHC. However, such predictions have to be confirmed by explicit trans¬ 
port simulations of the finite temperature and chemical potential QGP, since the macroscopic observables (Raa and V 2 ) depend 
not only on the microscopic processes and transport coefficients, as studied in this work, but also on the time evolution of the 
QGR Such studies will be carried out in the near future using the Parton-Hardon-String Dynamics (PHSD) transport approach 
[40, 57] for heavy-ion collisions from = 5 GeV to 10 GeV [58]. 


VII. SUMMARY 

We have presented in this work on- and off-shell approaches to describe the microscopic interactions between a heavy quark 
and the QGP degrees of freedom in a partonic medium at finite temperature T and quark chemical potential flq. In our formalism 
both the perturbative and non-perturbative parts of QCD are involved (DpQCD approach) and an off-shell description of heavy- 
quark interactions in the QGP (lEHTL) at finite T and jdq. For each of these models, we have calculated the differential cross 
sections and confronted their implications on the usual mesoscopic observables (energy loss, drag and diffusion coefficients, 
longitudinal and transverse Q momentum fiuctuations, etc.) at finite T and jlq. 

The fundamental parameters describing the heavy-quark collisional scattering have been fixed within the DQPM at finite 
temperature and chemical potential. The formulation of the DQPM at finite flq is based on the hypothesis that the phase boundary 
is close to the line of constant energy density in the {T^jiq) plane which leads to the approximation (II. 5). On the other hand this 
approximation is in very good agreement with recent IQCD results [51] on the expansion coefficient w.r.t. Accordingly we 
infer that the equation of state of QCD from Tc{llq) to higher temperatures is roughly under control at least for small/moderate 
quark chemical potentials. 

Our study demonstrates that even if the infiuence of the finite width of the quasi-particles on heavy quark scattering is small 
on the off-shell Q scattering cross sections, a noticeable effect is seen in the off-shell transport coefficients (lEHTL model) as 
compared to the on-shell ones (DpQCD model). This is due to reduced kinematical thresholds in the off-shell cross sections and 
is independent of the variables {T, jlq). 

A medium in which the chemical potential is finite leads to a reduction of the qQ and gQ elastic cross section and consequently 
to a reduction of heavy-quark energy and momentum losses as compared to a zero flq medium. Nevertheless, we have concluded 
that longitudinal momentum transfers are important not only in a hot medium but also in a dense medium whereas the dense 
medium leads to less transverse fiuctuations in the heavy quark propagation. The relative large drag at low temperatures in 
DpQCD/IEHTL is due to the strong increase of the running coupling as {T, flq) (infrared enhancement) for temperatures close to 

We have observed a smooth dependence of the energy loss on both variables T and jlq at finite but not too large values of 
jlq. Such a profile is expected for a cross-over transition from the partonic to the hadronic medium. For jlq = 0 the gluon mass 
depends on the temperature and therefore the increase of the energy loss is due to a change of the coupling. For jlq = 0.2 GeV, 
the energy loss is also increasing with temperature but less than for llq = 0 because here both the coupling and the effective gluon 
mass decrease and the increase of the infrared regulator is counterbalanced by the decrease of the coupling. Since the variations 
of all transport coefficients with T and }lq are rather smooth (within the present DQPM propagators) the transition from hadronic 
degrees of freedom to partonic ones remains a crossover up to jlq = 0.2 GeV. Present studies within the DQPM indicate a change 
to a first order transition only for > 0.3 GeV but this is yet model dependent and not robust. 

From the variations of the longitudinal and transverse momentum losses with T and jlq we expect a large suppression of Raa 
of charmed mesons in a hot and dense medium, but a much smaller value of the elliptic fiow V 2 in a hot and dense medium as 
compared to our knowledge from RHIC and LHC for finite temperature and approximately zero chemical potential. However, 
such expectations have to be confirmed by microscopic transport simulations (e.g. within PHSD) for heavy-ion collisions from 
AGS to SPS energies. This is also the energy regime of the future FAIR and NIC A facilities and the BES program at RHIC that 
all address the properties of QCD at high baryon densities or high flq, respectively. 
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